Using the X-ray Dust Scattering Halo of 4U 1624-490 
to determine distance and dust distributions 
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ABSTRACT 

: 

o ■ 

We present X-ray dust scattering halo results based on our 76 ks Chandra 
ACIS-S/HETGS observation of the LMXB dipping source 4U 1624-490. Through 
analysis of the halo light curves with 2-6 keV spectra over the persistent and 
\ dipping periods, we estimate a geometric distance of ~15 kpc to 4U 1624-490. 

We also fit halo radial profiles with different ISM dust grain models to assess the 
.^J \ location, uniformity, and density of the halo. Our analysis shows that the dust 

spatial distribution is not uniform along the line-of-sight; rather, it is consistent 
with the spiral arm structure mapped in Hn. The large difference between the 
absorption Hydrogen column (N^ s ~ 8 x 10 22 cm -2 ; probes all gas along the line- 
of-sight) derived from broadband spectral fitting, and the scattering Hydrogen 
column (Ajj ca ~ 4 x 10 22 cm -2 ; probes only Galactic gas) derived from our studies 
of the 4U 1624-490 X-ray halo suggests that a large fraction of the column is local 
to the X-ray binary. We also present (and apply) a new method for assessing the 
Chandra point spread function at large (> 50") angles, through use of the time 
delays from the observed dips. 



Subject headings: dust — scattering — X-rays: ISM — sources: 4U 1624-490 
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Introduction 



X-rays are not only absorbed but are also scattered by dust grains when they travel 
through the interstellar medium (ISM). This scattering of X-rays from a source behind a 
dust containing cloud in the ISM will lead to the formation of an X-ray scattering halo 
surrounding the X-ray source. The properties of the halo depend upon the size distribution 
and density of the dust grains, and on the relative geometry of dust, X-ray source, and 
observer. 



Such dust scattering halos were first d iscussed by lOverbeckl (119651 ). but they were not 
observationally confirmed until iRolJ (119831 ) observed the X-ray binary GX 339-4 with the 
Einstein X-ray Observatory. Since then, X-ray halos have been studied via observations 
facilitated by X-ray satellites which include Einstein, ROSAT, BeppoSAX, Chandra and 
X MM-Newton. Thus far, the most complete samples of these studies have been presented 
by IPredehl fc Schmittl Jl995h based on ROSAT data, and by Ixiang. Zhang fc Yaol (ho05h 
based on observations with the Chandra ACIS-S/HETGS. Of additional note. lDraine fc Tan 
( 120031 ) have analyzed ROSAT observations of the halo ass ociated with the X-ray nova 



VI 974 Cygni 1992 to confirm that the interstellar dust model of lWeingartner fc Drainei (12001 



hereafter WD01) is consistent with the observed X-ray halos. Further more, based on Chandra 

ACIS- I observations of the low-mass X-ray binary (LMXB) GX 13+l. lSmith. Edgar fc Shaferl 

(2002) reported that both the grain model WD01 and the "classical" model of lMathis. Rumpl fc Nordsieck 
( 119771 . hereafter MRN) can adequately reproduce its halo radial profile. 



X-ray halos are not limite d to the aforementioned b inaries. Most bright X-ray sources 
are surrounded by X-ray halos (IPredehl fc Schmittlll995l ). including one of th e most extreme 



of th e dipping sources, the Big Dipper 4U 1624-490, whi ch is an atoll source (ILommen et al. 



20051 ) with an IR counterpart of magnitude K s = 18.3 (IWachter et al.l 120051 ). The presence 
of a halo in 4U 1624-490 has been surm ised via BeppoSAX observatio ns, which exhibit a 
soft excess above its several arcmin PSF (jBaruciriska-Church et al.ll20001 ). The poor angular 
resolution of BeppoSAX, however, prevents the direct extraction of the halo radial profile at 
radii less than 100 arcsec, where the halo is much more prominent. 

In this paper, we present a focused study of the X-ray halo associated with 4U 1624- 
490, using the highest resolution (angular and energy) data to date, as afforded by the 



Chandra ACIS-S/HETGS. The ~3 hours long duration of the dips (jWatson et al. 



Smale. Church Balucihska- Church 2001; Balucihska-Church et al. 



1985 



200 ll ) and large obscu- 



ration of this compact object offer us the unique opportunity to use the time delay of photons 
arriving from the halo to determine the distance to 4U 1624-490 ($2] for theory and §4.21 for 
data analysis), and compare with dust grain models to assess composition and density along 
the line-of-sight (LOS, §4.3|) . (The large obscuration may be due to the accretion disk 
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stre am impacting the disk. A ~ 75% ob scuration has been reported by lWatson et al.lll985 
and I Church &: Balucihska- Church! Il995l ; the superior angular resolution of Chandra allows 
us to detect a 90% obscuration; §4.21 ) We also present a new method for assessing the 
Chandra Point Spread Function (PSF) at large angles ( <; 50"), and compare the PSF 
measured by our technique with Chandra Ray Tracing (ChaRT) predictions between 
9-160" ( EH) . 



2. Theoretical and Historical Background 

The theoretical calculations governing the observed halo surface brightness and the 
time delay of a s cattered photon with respect to an unscattered one have been discussed 



exten sively (e.g., iMauche fc Gorensteinl Il984l ; iMathis fc Led Il99ll ; iTriimper fc Schonfelder 



19731 ). Here, we briefly describe the main points. 



As discussed by lMauche fc Gorensteinl ( 11984 ) . the differential cross section in the Rayleigh- 
Gans approximation, coupled with the Gaussian approximation for a spherical particle of 
radius a, can be described by 



S(a, E, 9 sca ) 



d<J sca (a, E, 9 S 
dTl 



1.1 x 10 
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2Z 
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where K = 0A575(E/keV) 2 (a/ fim) 2 , and 9 sca is the angle of scattering. The mean atomic 
charge (Z), molecular w eight (M, in amu), mass density (p) and atomic scattering factor 
[F(E), from Henkel 1981 ] also factor into the calculation. A photon of wavelength A will be 
scattered to a typical angle 

9 = A/tt a oc E-\ (2) 

For small angles (several arcmin), we can approximate the observed angle by #~ (1 — x)9 sca , 
where x = d/D is the relative distance defined to be the ratio of the distance from scattering 
grain to observer (d) over the distance between source and observer (D); see Fig. [TJ 
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Fig. 1. — Geometry of the X-ray-scattering process for single (left) and double scattering 
(right). 
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The intensity of the observed first-order scattering halo, IscL(8,E), depends on the X- 
ray flux Fx{E) of the point source (4U 1624-490 in our case), and is a function of both 
the observed angle 9 and the energy E. The form of this equation, as initially derived by 



Mathis fc Led (119911 ). is expressed as 



li 1 c l(9,E)=F x (E)N n 



da n(a) / dxf(x)(l — x) x S I a, E 1 



1 — x 



(3) 



where iVn is the total hydrogen column density between the observer and the X-ray source, 
n(a) is the size distribution of the dust grains, and f(x) is the relative hydrogen density to 
the average total hydrogen density along the LOS at xD. For uniformly distributed dust, 
f(x) = 1, and eq. ([3]) can take on the more explicit form 



lM(e t E)<xF x (E)Nu 



da n(a) 



/ a V erfc(R) 
\fim J R 



(4) 



where R = K6 sca = K(6/l - x) and erfc(fl) = (2/^) jT dt exp (~t 2 ). 



Mathis fc Led (119911 ) further established that multiple scattering is likely important if 



there is enough scat tering optical dep t h to see an appreciable halo. As discussed by those 
authors, as well as by lPredehl fc Klosd (119961 ). for a fixed energy E doubly scattered radiation 
at the position of the observer can be described via 



li 2 c l(9,E) = F x (E) N u / dxf(x 



1 dx' f(x') 



\ — x 
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(5) 

(6) 
(7) 



For a typical scattering optical depth of r sca ps 0.5, the doubly scattered radiation 
dominates the multiple scattering terms at several arcmin such that higher order (> 2) 
scatterings can be largely neglected. 

Other factors, in particular the ISM dust composition and grain size distribution, affect 
our overall determination of the halo properties. Several different models exist for describing 
the composition and siz e distribution of ISM dust gr ains. The two most commonly used 
are the grain model of IWeingartner fc Draind (120011 . WD01), and the "classical" one of 
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Mathis. Rumpl &: Nordsieckl (119771 . MRN). The details of these two models are different, 



although they are both based on IR observations. The MRN model assumes both graphite 
and silicate grains with size distributions: n(a) oc a -3 5 , for a min < a < a max , where a min = 
0.005 /xm is the same for both grain types, whilst a max = 0.25 /xm for silicate grains and 
a max ~ 0.25 — 1 /xm for graphite grains. In contrast, the WD01 model, which includes very 
small carbonaceous grains (a < 0.005 /xm) and larger grains (a > 1 /xm), is comparatively 
more complex than the MRN model. To illustrate, the size of the carbonaceous grains 
extends to more than 1 /xm, while the number of grains decreases sharply with size - see 



Weingartner fc Draind (120011 ) for details. 



3. Observation 



We observed 4U 1624-490 on 2004 June 4 (MJD: 53160.26813, ObsID: 4559) with the 
Chandra High Energy Transmission Grating Spectrometer (HETGS) for 76 ks, covering one 
binary orbit. To reduce pileup, the observation was performed using a reduced 1.7 sec frame- 
time and one-half sub-array corresponding to 512 columns per CCD. Fig H] shows that our 
observation encompasses ~2.7 hr total of dipping periods (3 dipping events with durations 
respectively of about 3.5 ks, 2.3 ks and 4.0 ks). 



4. Data Analysis 



We used CIAO 3.3 with caldb 3.2 to extract HETGS spectra of the source for the 
persistent (non-dipping) and dipping periods (Fig.Hj). The 2—6 keV halo light curve between 
3" and 20" is extracted from the CCD S3 th order data, while the point source light curve in 
the same energy band is extracted from the dispersed data of the gratings (I s * order data). 
We note that the halo light curve, between 3" and 20", will have ~ 40% contamination 
from the bright point source photons scattered into the wings of the PSF. CIAO 3.3 is also 
used to extract the surface brightness of the th order data for the persistent and dipping 
periods. (Note that in order to ensure that neither the halo nor the PSF radial profiles 
are contaminated by the dispersed spectra of the HETGS, throughout we have excluded pie 
slices of the regions from the zeroth-order image that contain the MEG and HEG arms.) 
Our analysis technique for the results presented in §4.21 and §4.31 are best illustrated by the 
flowcharts of Figs. [3] and [5], respectively. 
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4.1. The broadband spectrum as determined from 1 st order HETGS data 

In order to determine the distance ( §4.2p and distributed Hydrogen column ( §4.31) along 
the LOS to 4U 1624-490, it is necessary that its light curve, and the total flux covering the 
entire ~ 76 ks spectrum, the dipping phase (~ 12 ks), and persistent phases (persistent 1 
+ persistent 2, ~ 59 ks) are estimated over the same energy bands. (This corresponds to 
2—6 keV over 200 eV incremental steps for this paper; see §4.21 for details.). Furthermore, 
because there is a direct proportionality between the source spectrum and the X-ray halo 
intensity such that any uncertainties associated with the spectral analysis will map onto the 
halo analysis, we take care to use spectra which suffer from the least amount of pile-up in 
our continuum modeling. 



4- 1.1. The time averaged continuum covering the binary cycle 

The best fit model based on HEG±1 fits to the time-averaged ~ 76 ks spectrum between 
1.5-10.0 keV is 1.371^5 keV blackbody plus V = 1.0±g:| power-law modified by JV H = 

7 - 6 -o:i x 1022 

cm . Not surprisingly, these parameters are similar to that derived from the 
persistent period which is discussed in more depth subsequently. 



4-1.2. The persistent phase continuum 

The maximal count rate for the HEG±1 spectra during this phase is 7.8 x 10~ 3 counts 
per pixel per frame time, corresponding to a maximal pileup ~ 2% at ~5 keV. The maximal 
count rate per pixel in the MEG±1 is about two times that of HEG±1, indicating a maximal 
pileup of ~ 4%. Therefore, we have taken care to use the nearly pileup free HEG±1 spectra 
to derive the flux for the different bands used in our analysis of §4.21 and §4.31 The halo 
intensity is proportional to E~ 2 such that it decreases sharply with increasing energy. We 
therefore retain as many low energy (~ 2 keV) halo photons as possible for our halo studies. 
We further note that we derive broad band spectral parameters based on the 1.5-10 keV 
spectral region, whilst only the 2-6 keV flux (measured in 200 eV steps) is used in order to 
match our halo studies of §4.31 

Using0xsPEC 12.2.1 ( Arnaud 1996 ). we fit the broadband source continuum HEG±1 



spectrum of the persistent phase with various combinations of powerlaw, blackbody, diskbb 



1 http://heasarc. gsfc.nasa.gov/docs/xanadu/xspec/index. html 
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and thermal bremsstrahlung m odified by the Tuebingen-Boulder ISM absorption model 
( IWilms. Allen &: McCrayl |2000| ). (Because of the high count rate of the source spectrum, 
the data were binned to require at least 400 counts in each energy bin.) The fitting results 
are listed in Table [TJ We find that both Model 1: absorption x (blackbody+powerlaw) and 
Model 2: absorption x (diskbb+power law) fit the data well. However, the latter gives an un- 
reasonabl e photon index of T = 12j j j, whi lst the former agrees better wi th the parameters re- 

ported by lBalucihska-Church et al.l (120001) based on Beppo SAX data, by lSmale. Church fc Balucihska-Chun 
(120011 ) based on RXTE data, and by lParmar et al.l (120021 ) based on XMM-Newton data, bar- 
ring a flatter photon index in our fitting. The steepness of T appears strangely tied to choice 
of cold absorption model (i.e. WABS versus tbabs), which we will explore in more depth 
in our forthcoming paper on the high resolution Chandra spectrum of 4U 1624-490. For 
present purposes, it suffices that Model 1 describes the broadband continuum spectrum well 
(see Fig. [2]), such that we can confidently extract flux values over the incremental 200 eV 
steps between 2-6 keV for our calculations of §4.31 



4-1.3. The dip-phase continuum 

We also fit the MEG±1 and HEG±1 spectrum of the dip phase with many of the same 
models. Due to the lower overall count rate during this phase, we were only able to rebin 
these data to require S/N~8 bin -1 . Both the MEG±1 and HEG±1 dip spectra are pileup 
free. We find the best fit for this broad band dip continuum spectrum to be a T = 0.7 ± 0.5 
power- law modified by JV H = 9.05^ x 10 cm cold absorption. Fluxes in 200 eV steps 
were also obtained in the same 2—6 keV energy band to compare with the persistent periods. 



4.2. Using the halo and point source light curves to determine distances 



The scattering photon travels longer distances than the unscattered one. The delay 
time dt is given by 



dt = 1.15 h (-D/1 kpc) (9/1 arcmin) 



x 



1 — x ' 



Thus measurement of a time delay between the point source and the halo yields direct 
information about the distance to the dust cloud and the source. The ~ 3 hr long dipping of 
4U 1 624-49 provides an excellent opport unity to search for such a delay. iTriimper &: Schonfelder 
(119731 ) and Ku. McCray &: Kelleyl (119861 ) propos ed methods fo r using this behavior to mea- 
sure the distance to variable X-ray sources, and IPredehl et al.l (120001 ) presented a successful 
geometric distance determination to Cyg X-3 based on Chandra ACIS-S/HETGS observa- 
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Table 1. Spectral fits to non-dip continuum. 



Model 


Nn (10 22 cm" 2 ) 


kT (keV) 


r 


X 2 /dof 


blackbody + Power law 


7 s+ - 8 

'•S-0.6 


1 Q4+0.09 


i 1+0.5 
■ L - A — 1.0 


364/379 


Disk blackbody + Power law 




2.57±0.06 


12±1 


356/379 


Disk blackbody 


8.76±0.13 


2.64±0.05 




373/381 


Bremsstrahlung 


10.02±0.14 


10.318-J 




442/381 


Blackbody 


6.60±0.11 


1.54±0.02 




428/381 


Power law 


10.77±0.18 




1.90±0.04 


510/381 



Note. - These results come from fits to the Chandra ACIS-S/HETGS HEG±1. The 
MEG±1 is not used since it suffered from more pileup than HEG±1. Errors are 90% con- 
fidence. All models listed are m odified by the Tuebingen-Boulder ISM absorption model of 



Wilms. Allen fc McCrayl (120001 ) . as noted in the iV H values. 




Fig. 2. — The best fit models over-plotted on the broad band spectra of 4U 1624-490 during 
the persistent phase (HEG±1; left) and dipping periods (HEG±1 and MEG±1; right). 
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tions. 

The time delay dt is only dependent on the location of scattering, but the halo surface 
brightness is determined by the size, position, and composition of the dust grains, as well 
as the source flux and scattering hydrogen column. Therefore, the spectrum of the source, 
the dust grains model and the spatial distribution of the grains is implicitly factored into 
eq. ©. 

The halo intensity is proportional to E~ 2 such that it decreases dramatically at the 
very high energies. Fig. [2] shows that the majority of the 4U 1624-490 counts are contained 
within the 2-6 keV energy band. Therefore, we restrict our comparative analysis of the halo 
and point source light curve exclusively to this band, where fluxes are estimated over 200 eV 
steps based on the continuum model described in §4.1.11 

In order to avoid pileup effects, the bright point source light curve is extracted directly 
from the dispersed 1 st order data of the HETGS, where the count rate is a factor of ~ 4 lower 
than that of the non-dispersed th order data. In contrast, we extract the (lower flux) halo 
light curve between 3" and 20" from the th order data, where the minimum angle is restricted 
to 3" to mitigate pileup. We estimate that while pileup at 3" is less than 2%, between 3"-20" 
it is <C 1% so that we can safely neglect it for our analysis. The maximum angle restriction 
of 20" was chosen to reduce multiple scattering effects, which if not factored in properly can 
lead to over-predictions of the time delay. Alternatively, we could have restricted the energy 
band to 3-6 keV, for example. However, this not only would have reduced the degree of 
multiple scatterings, it would have greatly decreased the halo intensity, which scales as E~ 2 . 
Thus we opted for the angular restriction with an energy that covers well the observable point 
source continuum. We estimate that the total halo intensity between 3"-20" associated with 
the second-order scattering is less than 5% that of the first-order scattering in the 2-6 keV 
energy band of interest. While we account for 2 nd order scattering effects in our analysis, we 
wanted to make sure that it would not play a major role in affecting our final results. 

A comparison of the halo and point source light-curve (Fig. H|) reveals a mean delay 
of ~1.6 ks for source photons arriving from the halo that is at 3"-20" from 4U 1624-490. 
It is interesting to note that the dips we detect show as much as ^ 90% blockage of the 
compact object (Fig . HI), com pared, for example, to t he 75% report ed previously in EXOSAT 



(IWatson et al.lll985T). Ginaa (I Jones fc Watson! Il989l). BevpoSA X ( iBaluciriska- Church et al. 
200ol ) and RXTE ( Smale. Church fc Balucihska- Church 2001 ) observations. This is at 



tributed to the fact that the superior imaging capabilities of Chandra can better separate 
the source light-curve from contaminated light from the halo. Using the time delay, and 
eq. (jU), a simple estimate for the distance to 4U 1624-490 can be derived. If we assume a 
mean value for the fractional distance: x — | and a mean effective angle 6 = 10.3" derived 
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Chandra HETG/ACIS-S 
Data of 4U 1624-490 



1 st Order Data 



Flux in different energy 
bands from HEG ±1 

f t0 J E > 



A 

Point Source 
Light Curves 
LC (t) 

pomV 



Assumed 
Distance 
D 



Light Curves with the 
Dust Grains Model 
(WD01 and MRN) 



th Order Data 
(in CCD S3) 



Spectra HEG ±1 


3 






and MEG ±1 




5 




Total S(E) 







f 

Halo 
Light Curves 

halo 
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Getting Minimalx 2 



Distance of Source 



Fig. 3. — Flowchart showing the data analysis process for determining the distance to 
4U 1624-490. 
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Fig. 4. — The 2-6 keV halo and point source light curves. The top-most dotted line is the 
point source light curve for 4U 1624-490. The combined halo and PSF light curve between 3" 
and 20" is shown, with error bars, in the middle. The overplotted solid line is the "best fit" 
(x 2 statistics) for the halo light curve based on LC WD01 (t) (dashed line) convolved with the 
PSF lightcurve. To facilitate clearer viewing, we multiply the observed halo plus PSF light 
curve by a factor 4, the LC WDm (t) from model WD01 by a factor of 3, and the LC M RN(t) 
from model MRN by a factor of 2. 
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using 



9d6 exp 



-Ka 2 E 2 



1 — x 



exp 



-Ka 2 E 2 



1 — x 



(9) 



where 9 m i n = 3 arcsec and 6* ma3; =20 arcsec, we obtain a distance to 4U 1624-490 to be about 
13.1 kpc. In the subsequent discussion, we discuss a more rigorous approach for deriving the 
distance. 



Following the prescription presented by lPredehl et al.l (120001 ). we initially generate model 
light curves LCwdoi(^) (using the WD01 dust model) and LCmrn(^) (MRN dust model) for 
the halo (Fig.HJ), according to eqs. (jSj) & ((HJ), assuming a uniform spatial distribution of dust, 
for an initial guess of the distance D=10 kpc. [LCwdoi(^) an d -^Cmrn(^) will hereafter be 
referred to generically as LC mo dei(t).} Then we use LC mo dei(t) together with the light curve 
of the point source LC po i nt (t) to fit the light curve of the halo LChaioit), allowing for the 
ratio, LC mo dei (t):LC p i n t (t) to vary until a minimal \ 2 is found. This is iterated for input 
distances ranging between 8—20 kpc in 0.1 kpc steps. The best fit distance to 4U 1624-490 
derived from the WD01 and MRN models for uniformly distributed dust are respectively 
-Dwdoi = 15. kP c and -Dmrn = 10. 2^ 1 ; 4 kpc. (Errors are quoted at 90% confidence 
based on variations in x 2 with D.) Within errors, both models appear to give consistent 
values for the distance. However, it may be of interest to note that the distance derived 
from the MRN model is somewhat sensitive to the allowed maximum size a max for the dust 
grains. Since the maximum size for carbonaceous grains in the MRN model is not specified, 
we tried to test our distance determination by changing the maximum size for carbonaceous 



grains in the MRN model. For a max = 0.25 /im, we obtain -Dmrn = 10.2_ 1 ; 4 kpc, as noted 
above, in comparison to -D M rn — 15.0jJ;g kpc for a max = 0.42 /im grain sizes. The latter 
distance estimate is comparable to the value obtained based on the WD01 model. 



The W DOl model likely provides a m ore realistic present day picture of the dust grains 
in the ISM. ISmith. Edgar fc Shaferl (120021 ). however, reported that the MRN model is better 
than the WDOl model in fitting the halo radial profile of GX 13+1, based on observations 
using the Chandra ACIS-I. R. Smith (priv. comm., Dec. 2006) notes, however, that this 
may be attributable to systematic errors. In a recent reanalysis of GX 13+1, based on a 
new Chandra calibration, Smith finds that the WDOl model may be preferred (although 
the differences between the WDOl and MRN models are still small). Therefore, for our 
distance determinations, we are still inclined to accept the value of -D 4U i624 = 15.0jl|;| kpc 
ba sed on fits using the WDO l model. This is consistent with the 10-20 kpc value derived 
by IChristian fc Swank! (119971 ). based on a method which compared the hydrogen column 
densities from their spectral fitting to an exponential distribution model of hydrogen in the 
Galaxy. 
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Furthermore, since dust grains between us and the point source are unlikely to be uni- 
formly distributed, we check whether and to what extent our distance determinations are 
affected by unevenly distributed dust. Accordingly, in our modeling, we vary the dust place- 
ment in additional fits, e.g., x=0.5-l, 0.4-1, etc., where x=0 and x=l correspond respectively 
to our position and that of the source (see Fig. [8]). We find our distance estimates to be robust 
to these changes, as long as the dust is uniformly distributed near the source. This is because 
the distance estimates are based primarily on small angle scatterings , where the observe d 
halo is attributed primarily to dust near the source (see, e.g., Fig 4 of lMathis fc Ledll99ll ). 
This finding is consistent with our results based on fits to the halo profile for determining 
column densities in §4.31 



4.3. Using the halo radial profile to determine the A^ H spatial distribution 

Fits to the halo light curve show that ~40% of the photons come from the point source 
as far out as 20". While Chandra has very good imaging resolution, the halo brightness is 
not much greater than the PSF (see Fig. right). Therefore, a good understanding of the 
telescope PSF is necessary for optimal results. 



4.3.1. The Chandra PSF and halo of 4U 1 624 -4 90 



As discussed by ISmith. Edgar fc Shaferl (12002 ). the CIAO tool mkpsf and the Chandra 



raytrace model SAOsac recreate the observed core of the PSF well but underestimate the 
wings of the PSF. The same is true of the tool ChaRT since it also uses SAOsac. Column 
(d) of Table [2] lists the comparison between the PSF (PSFchaRr) from the ChaRT simu- 
lation and the one (PSF HerX i) derived from the Chandra AC1S-S/HETGS observation of 
the almost halo-free source, Her X-l (ObsID 6149). Parameters derived from the Her X-l 
observation are used in its ChaRT simulation. As such, we can confidently use ChaRT to 
assess the PSF behavior at small ( ^ 50") angles. For angles larger than 50", we subse- 
quently present a new method for assessing the Chandra PSF by using the observed dips in 
4U 1624-490. 

Fig. [6]shows the WD01 halo model light curve over different angular distances compared 
against the point source light curve. It is clear that at increasingly larger angles, the flux 
behavior of the point source is reflected less in the halo light curve (Table [2]). This is in 
part due to a "smearing" effect at large angles as well as the increased delay between the 
photons arriving from the point source and those arriving from the halo. Capitalizing on this 
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Fig. 5. — Flowchart showing the data analysis process for determining the PSF radial profile 
and halo radial profile which are used to assess distributions of LOS N H . 
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2.0273 2.0274 2.0275 2.0276 2.0277 2.0278 

Time (10 8 s) 

Fig. 6. — The theoretical halo light curves derived from the WD01 model (dashed) at different 
angular distances based on the observed point source (4U 1624-490) light curve (solid). The 
halo light curves are arbitrarily multiplied by some constant factor for easier viewing. 
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observed effect of the dipping phenomenon on the halo light curve, we discuss a prescription 
for assessing the Chandra PSF for 4U 1624-490 at large angles. First we define: 

PSF(r) = [I persistent (r) - I dip {r)\ F "*«*** (10) 

\" persistent "dip) 

where Idi p (r) and I per sistent( r ) are respectively, the surface brightness over different radii 
during the time periods when 4U 1624-490 is dipping and when it is not. Similarly, i^ ip and 
Fpersistent refer to the point source flux during these period^. Based on this, we can then 
define the surface brightness of the halo : 



-^halo(^) I per sistentir*) PSF{T) \1 dip(,T^) ^persistent Ipersistenty^) F d ip\ / '[Fpersistent Fdip)- 

(11) 



To track properly the spectra of the dip and persistent periods and determine how these 
spectra are reflected in the behavior of the energy dependent PSF, we use eq. ffTUl) to calculate 
the PSF over the energy range 2-6 keV in 200 eV steps. Note that this equation implicitly 
assumes no time-dependence of the halo flux. This is, of course, inaccurate, especially at 
small angular distances from the source. It becomes more accurate, however, at large angular 
distances where the dipping behavior is "smeared out" in the halo lightcurves (see Fig. 
Thus, we expect eq. (1T01 to be less dominated by systematic uncertainties at larger angular 
distances. Large angular distances are also where ChaRT estimates of the PSF become 
more problematic. 

To assess the systematic errors in our determination of the PSF radial profile, we use 
eqs. (J3J) and (jSJ) to calculate the theoretical halo light curves [LC mo d(t)] at different ranges 
of angular distances, as shown in Fig. [61 Specifically, as described in §4.21 we used the 
lightcurve of the persistent source with different presumed models for the source distance 
and dust distribution, and chose a model that minimizes the x 2 f° r the fit to the halo 
lightcurves. The halo intensity averaged over the dip phases [i.e., essentially the estimate 
embodied in eq. (11 ID ] can than be compared to the halo intensity in the persistent phases 
(i.e., the "true" steady state halo flux). This estimate is presented in column (a) of Table 
[2j As expected, for radii near the point source, selecting the periods during the dips slightly 
underestimates the persistent halo flux, and hence leads to a systematic overestimation at 
small angular radii of the PSF as calculated via eq. (JTD]) [column (c) of Table [2J. 



2 As the halo dips affect the halo surface brightness, especially shortly after the dips, throughout we 
calculate I per sistent solely from persistent phase 1, i.e., the w 20ksec before the dips. However, for the point 
source flux, F pers i ate nt, we use both persistent phases 1 and 2. 
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As shown in Table [21 for radii < 50", ChaRT yields less than 13% uncertainties for 
the PSF, while for radii > 50", using eq. ffTUl) yields uncertainties < 10% for the PSF. The 
relative uncertainty in the halo flux is approximately given by the fractional uncertainty in 
the PSF multiplied by the ratio of the PSF flux to the halo flux [i.e., columns (c) multiplied 
by column (b) or columns (d) multiplied by column (e) from Table [2]. Thus, by using 
the ChaRT PSF estimate for radii < 50", and eq. (jTUl) for radii > 50", the systematic 
uncertainty in the halo profile should be < 3% everywhere. 



4-3.2. Fits to the halo radial profile 



Having gained some understanding of the PSF behavior, we us e the WD01 and MRN 



model s to fit the halo radial profile using the fitting codes develo ped bvlSmith. Edgar fc Shaferi 
( 2002 ). as also applied to the Chandra halo studies presented in Xiang. Zhang fc Yaol ( 2005 ). 



See also Fig. [5] for a flow chart of the analysis process. We confine our halo fitting to angular 
distributions corresponding to radii 9" to 160" where the extreme of pileup effect is < 0.5% 
at ^ 9" and decreases with angular distance. Therefore, we can assume a pileup free halo 
for our analysis. 

In order to assess dust properties and distribution, for our initial fitting, we assume 
unifo rmly distributed dust be tween x=0 (us) and x=l (4U 1624-490; see Fig. [HI as recreated 
from lCaswell fc Hayneslll987l and annotated to show the distance of 4U 1624-490 from us). 
Based on this scenario, we find that neither the MRN model (with a max (graphite) = 0.42 /zm; 
x 2 /dof = 141/49 for a best fit N H = 4.44+°;^ x 10 22 cm" 2 ) nor the WD01 model (x 2 /dof = 




Fig. 7.— PSF (left) and halo (right) radial profile of 4U 1624-490. 
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Table 2. Estimation of the Systematic Error of PSF and Halo. 



radius (") 


jwd01,dip, ) ™ d 01( r ) a 
halo v ' halo V > 


PSF data (r) b 


APSF data (r) c 


PSF HerX1 (r)-PSF ChaRT (r) d 


PSF ChaRT (r) e 




T Zio( r ) 


PSF data (r) 


PSF HerX i(r) 


'halo VJ 


09-20 


-0.12 


0.77 


0.18 


0.02 


0.36 


20-30 


-0.07 


0.51 


0.16 


0.08 


0.23 


30-40 


-0.05 


0.46 


0.13 


0.13 


0.20 


40-50 


-0.04 


0.39 


0.11 


0.13 


0.19 


50-60 


-0.03 


0.32 


0.09 


0.24 


0.16 


60-70 


-0.02 


0.24 


0.09 


0.30 


0.13 


70-80 


-0.01 


0.36 


0.04 


0.32 


0.11 


80-90 


-0.01 


0.25 


0.04 


0.38 


0.11 


90-100 


-0.003 


0.16 


0.02 


0.46 


<0.10 


100-160 


< -0.003 


-0.30 


<0.01 


> 0.40 


<0.10 



a J^° 1 (r) and I™ alo ' %p {r) mean the halo intensity for the persistent and (averaged) dip 
phases, respectively, derived from fitting the energy- dependent lightcurves with the WD01 
model. 

h P S F data {r) is the PSF intensity derived from eq. (IIP]) , while I^„(r) is the halo intensity 
derived from eq. (TTTT) . 

c APSF data (r) is an estimate of systematic error in the (persistent phase) PSF due to the 
underestimate of the derived halo intensity [i.e., column (a)]. Specifically, we use the estimate 
APSF data (r)/PSF data (r) « [[Iffi(r) - lZ^(r))/Iffl(r)) [F persistent /(F persistent - F dip )} 
[ltlo(r)/PSF data (r)]. 

d (PSF HerX i — PSFchaRr) I ' PSF HerX i means the systematic error of PSF intensity from 
ChaRT simulations compared with the Her X-l observational data. 

e Ihaio RT ( r ) = Ipersistentij) — P S FchaRT{r) is the halo intensity with the ChaRT simualtion 
PSF. 

Note. — The data in this table is is used solely to discern the relative accuracy of our PSF 
estimates, and are not used in the calculation of either the PSF or halo radial profiles. 
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171/49 with JV H = 3.40i°;$j x 1q22 ™~ 2 ) provides good fits to the halo profile. Both models 
under-predict the halo surface brightness at >80" by ~50%. We also fit the halo radial 
profile to include the second-order scattering and find a 50% substantial improvement in our 
fits. (We remind the reader that the angular regions probed in this study are much larger 
than the distance determination study of §4.21 where multiple scattering effects were found 
to be negligible.) While the fits improved somewhat, they are still statistically unacceptable 
with x 2 /dof ~ 91/49 for the WD01 model, and x*/dof ~ 89/49 for the MRN model of 
a max (graphite) = 0.42 /zm. 

Having determined that the uniformly distributed dust between us and the point source 
4U 1624-490 is an unlikely scenario, we investigate whether the models are sensitive to 
patchy distributions. Accordingly, based on Fig. [HI we roughly divide the dust distribution 
along our LOS into three parts : x = 0.0-0.20 (Region 1: Rl), 0.20-0.40 (R2), 0.40-1.0 (R3), 
corresponding to a distance d relative to us of, respectively, 0—3.0 kpc (Rl), 3.0—6.0 kpc 
(R2), and 6.0—15.0 kpc (R3), as seen in Fig. [9j While dust is assumed to be smoothly 
distributed in each of these regions, the quantity is allowed to vary independently. Based on 
this, fits to the halo radial profile using the WD01 model yield Hydrogen column densities for 
each region (Table[3J. Compared to fits assuming uniformly distributed dust between x=0— 1, 
Ax 2 = 26 for 47 d.o.f. (xl — 1-38), i.e., >99% confidence for 2 additional parameters, 
according to the F-test. Similarly, fits based on three parts with the MRN model give 
Ax 2 = 24 for 47 d.o.f. (xl — 1-38), i.e., again >99% confidence for 2 additional parameters. 
The second-order scattering is accounted for in these fits by adding the numerically integrated 
value of the second-order scattering to the first-order scattering value during fitting of the 
halo radial profile. The scattering hydrogen column density #h mrn = 4.8 x 10 22 cm -2 derived 
from the MRN model i s consistent with the resu l t (Nw = 5.0 x 10 22 cm~ 2 corresponding to 
t = 2.4) obtained by iBaluciriska- Church et al.l (120001 ). We also notice, however, that a 



smaller r = 1.8 obtained by iTrigo et all (120061 ) is more consistent with the result A^ HwD01 



3.6 x 10 22 cm 2 derived from the WD01 model. 

The derived hydrogen column density ^Hwdoi = (1-58±0.15) x 10 22 cm -2 (translated into 
a particle density n^i ~1.7 cm~ 3 ) for the two nearby spiral arms enco mpassed within Rl, is 



consistent with the result derived from the X-ray halo of Circinus XI by lXiang. Zhang fc Yao 



Table 3. Halo fits for deriving Nn based on a three parts dust distribution scenario. 



Model 


A^IO 22 cm" 2 ) 


JV« 2 (10 22 cm" 2 ) 


AT« 3 (10 22 cm" 2 ) 


^'(lO 22 cm" 2 ) 


X 2 /d.o.f 


WD01 
MRN 


1.58 ±0.15 
9 n«+ - 19 

^■ U °-0.20 


00+ 18 
u - uu -o.oo 

u - uu -o.oo 


2 04+ 08 
z - u<t -0 .09 

2.73 ± 0.11 


n 09+O.O6 
°- oz -0.06 
, s1 +0.08 


65/47 
65/47 
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Fig. 8. — The spiral arm structur e of the Milky Way Galaxy and the position of 4U 1624-490 
in it. The figure is recreated from lCaswell fc Hayned (119871 ). The positioning and the source 
is based on the distance derived in §4.21 
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Fig. 9. — X-ray halos observed at 2.0-6.0 keV, fit with dust grains model WD01 (top) and 
MRN (down). The solid line is the model where the dust spatial distribution is divided into 
three parts, which is marked in the Fig. [8j and the dotted line is the model with uniformly 
distributed dust, respectively. The dashed line is the double scattering intensity model with 
3 parts dusts. 
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( 2005 ) over similar distances (A'h^ 1 = 1.4 x 10 22 cm -2 at ~ 3 kpc from us) along a similar 
LOS (/ = 322°, b = 0.04°). It should be noted that these two values do not need to be exactly 
equal since the LOS to 4U 1624-490 and the one to Circinus X-l pass through different regions 
of the same two spiral arms. For R2 (the region between spiral arms 2 and 3), we estimate 
an upper limit to the particle density ~0.2 cm -3 ) based on a 90% confidence value for 
iVjj < 1.8 x 10 21 cm -2 (Table EJ). Combined, the hydrogen particle density in Rl is nju ~l-7 
cm -3 , twice as high as the average density along the entire LOS to 4U 1624-490. The best 
halo radial profile fits are shown in Fig. [9J 

While the absorption hydrogen column density includes both the Galactic ISM and the 
gas nearby the source, the scattering hydrogen column density only comes from the Galactic 
interstellar medium. Therefore, our finding that the scattering A^h derived independently 
from the WD01 and MRN models is much less than the absorption Nr derived from spectral 
fitting, implies that there is significant absorption intrinsic to the source. If we naively take 
the difference between the best fit value N^ s ~ 8 x 10 22 cm -2 of the broadband HETGS 
spectra due to Galactic and source ISM and N^ a ~ 4 x 10 22 cm -2 derived from our halo 
studies, a Ah ^ 4 x 10 22 cm~ 2 can conceivably be local to the source, possibly put there 
by the stellar wind of the companion. 



5. Summary 

We improve upon previous distance estimates for 4U 1624-490 by making use of the 
delay time of the halo photons relative to the bright point source photons to obtain 
-D4 TI1624 = 15.0^ 2 fi kpc . This is consistent within the errors of the 10-20 kpc estimates 



by IChristian fc Swank! (119971 ) using a different technique 



We find that varying dust distributions will not affect our distance determination to 
4U 1624-490 except possibly for a scenario where there is no dust within ~7.5 kpc of 
the source (x=0. 5-1.0). This extreme scenario does not match the Ah derived from 
our halo fitting results so it can be ignored as potentially problematic for our distance 
estimates. 

Using the extreme dipping behavior of 4U 1624-490, we discuss a new method for 
estimating the Chandra PSF at large angles (> 50"). In a comparison with ChaRT 
estimates for the PSF, we find that if we estimate the PSF using ChaRT at < 50", 
and our method at > 50", we can limit the errors associated with our halo analysis to 
& 3% over the angular 9"— 160" region. 

Varying dust distribution does affect the derived column densities. A simple estimate 
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based on our halo fits imply the hydrogen particle density in the spiral arms is tiri ~1.7 
cm -3 , and the one between two spiral arms <0.2 cm -3 . 

For the future, larger field-of-view and high throughput observations combined with 
on-going Chandra studies will allow us to better diagnose the scattering of dust near 
and far from us to reveal more detailed spatial information. 



Note : After submission of this paper, llaria et al.l (120061 ) posted a preprint presenting the 



spectral analysis of our HETGS data for 4U 1624-490 which is currently public. Our complete 
analysis of these data will be presented in a forthcoming publication. 

We wish to thank Randall Smith for insightful suggestions. This work was funded by the 
NASA / Chandra grant GO4-3056X - we are thankful for its support. We are also grateful 
to the Harvard University Clark fund for research support. 
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